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Abstract: Bloch branes were introduced previously and are constructed in a system 
described by two real scalar fields coupled with gravity in (4, 1) dimensions in warped 
spacetime involving one extra dimension. This work investigates gravity on such thick 
branes with internal structure and focuses on the effects of massive graviton modes 
localization on the brane and to what extent they might reproduce the 4d gravity at a 
scale before escaping into the extra dimension. In this way gravitational measurements on 
the brane could reveal the existence of the extra dimension on some scales, with possible 
applications on brane cosmology. 
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1. Introduction 

Domain walls as extended defects in field theory have been used in high energy physics 
for representing brane scenarios with extra dimensions |jl], ||, ^. Considered as domain 
walls with an internal structure, the Bloch brane is constructed in a (4, 1) model of 
two scalar fields coupled with gravity. The interaction among the fields depends on a real 
parameter, which generalizes the standard (/>^ model, changing the way the scalar field self- 
interacts. In this way the brane is generated dynamically and the asymptotic bulk metric 
is a slice of an AdS^ space. This and other thick brane |^, ^ models are alternatives to the 
infinitely thin brane models (see, for instance, Ref. |^ and extensions), there constructed 
with non-dynamical source terms in the action. 

In a previous work the internal structure of the Bloch brane was analyzed in terms 
of the matter energy density as a function of the extra dimension. It was found that in a 
certain region of the coupling parameter, an splitting of the defect in two interfaces occurs. 
In that work the effect of gravity on the internal structure of the domain wall was made 
evident as an induced attraction between the internal interfaces that form the defect. 

This intriguing possibility has led us to refine the former investigations and search for 
resonance states Q for our r small. We found that this occurs in the region of parameter 
where the brane splitting is more accentuated as we report on section 3. In the same 
section we considered the Newtonian potential for the whole range of the r parameter with 
the aim to study gravity localization. This is the main issue of the present work. Before 
studying this problem, however, let us first review some of the results of the Bloch brane 
model. It is described by the action 



S = I d^xdy^M[-\R+lda(t>d''^ + ^daXd''x-V{^,x) 
where g = det{gab) and the metric 



fl.l) 



ds"^ = dabdx^'dx^ 
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describes a background with 4-dimensional Poincare symmetry with y as the extra dimen- 
sion. Here a, 6 = 0, 1,2,3,4, and e^^ is the warp factor. We suppose that the scalar fields 
and the warp factor only depend on the extra coordinate y. 

The action given by Eq.(^T^) leads to the following coupled differential equations for 
the scalar fields (j){y)-, xiv) the function A{y) from the warp factor: 



(f)" + 4A'(t>' 
A" 



A 



12 



dx 



+ X'^)-Iv{ct>,x) 



(1.3) 
(1.4) 
(1.5) 
(1.6) 



where prime stands for derivative with respect to y. 
The potential used is 



mx) 



1 



dW 



+ 



dW 
dx 



where 



T^((/>,X) = 20--(/>3-2r(/.x^ 



(1.7) 



(1.8) 



and r is a real parameter. The particular relation from Eq.( |l.7| ) between V and W leads to 
a description of the scalar fields in terms of a coupled set of first-order differential equations, 
with analytical expressions for 0(2/) and xiu) identical to the obtained in Refs. |9| for 
flat spacetime. 

The first-order differential equations which also solve the equations of motion are 



1 dW_ 

2 "5^ 



(1.9) 



and 



1 dW 



A' = --W. (1.11) 

3 ^ ' 



Similar investigations where also done in Ref. [^]. 

In the present case the former equations can be solved analytically leading to the 
following results 

ct>{y) = tanh(2ry), (1.12) 
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X{y) = y - 2 ) sech(2r2/) (1.13) 
and 

Aiy) = —\{1 - 3r)tanh^(2ry) - 2 In cosh(2ry)j (1.14) 
9r L J 

Graphical analysis 0] from these solutions for several values of the parameter < r < 0.5 
show that the thickness of the defect increases with decreasing r, with A{Q) = implying 
canonical normalization for the 4-dimensional metric in the transverse space on y = ||l3| . 
Note also from Eq. ( [L.13 ) that for the limit r = 0.5 the one-field scenario is recovered. 



2. Fluctuations and massive modes 

We consider now the effective 4-dimensional gravitational fluctuations in the conformally 
flat background discussed previously, as well as that of scalar flelds around solutions 
(Eqs- CTra ): 

ds^ = e^^^y^ + ehf,^)dx''dx'' - dy'^ (2.1) 

and we set (j) ^ (j) + ecj) and X ^ X + ^X-, where h^j^y = h^y{x, y), (p = (/)(x, y) and X = x{x, y) 
represent small perturbations. On the transverse traceless gauge the metric perturbation 
separates from the scalars [|^], leading to 

hl, + AA'h'^, = e~^^nh^,, (2.2) 

with □ the (3, l)-dimensional d'Alembertian. 

This equation can be decoupled by separating the 4-dimensional plane wave pertur- 
bations in a form of plane waves from the extra dimension contribution. We introduce a 
new variable z that turns the metric into a conformal one. This changes the equation for 
the extra dimension contribution of the metric perturbations in a Schroedinger-like form, 
where no single derivative terms are present. The new conformal coordinate z is deflned 
by dz = e~^^y^dy. The separation of variables taken as 

V(x,z) = e*P-^e-t^(^)V/..(^), (2.3) 

turns Eq. (|2.2| ) into a Klein-Gordon equation for the 4-dimensional components of the 
transverse-traceless hfj_^, with the remaining Schroedinger-like equation 

- ^-^^^ + Vsch{z)Mz) = m^^miz) (2.4) 

for its fifth-dimensional component, where we dropped the fj.i' indices for the wavefunctions, 
now labeled by their corresponding energy m?. Here the potential is given by 

Vsck{z) = lA"iz) + lA'\z) (2.5) 
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Particularly interesting is that in this new conformal variable the appropriate inner-product 



for the wavefunctions ip{z) is the conventional quantum mechanical one |13|; in this way 
we can interpret a normalized \^jJmiz)\'^ as the probability for finding a massive mode at 
the position labeled by z. 

The solutions of Eq.(2.4) for general m form a tower of states associated with the 
massive modes responsible for the gravitational interaction in the semiclassical theory we 



are considering. Note that the l.h.s. of Eq.(2.4) can be written as a positive definite 
Hermitean operator, meaning absence of tachyonic excitations |14] with a zero mode (nor- 
malized (m = 0) state) given by 



(2.6) 



where iV is a normalization factor. 

Graphical analysis from Eq. (|2.5| ) shows that for some region of parameters Vsch{z) is 
a "volcano-type potential" with a true minimum at y = for r* < r < 0.5, whereas for 
r < r* there appears a splitting of the potential with two separate minima. In Fig. I we 
reproduced the figure from Ref. Q in order to ease the comparison with the results from 
next section. As the potential goes to zero as \z\ — > ±00, the possibility for the model 
to reproduce 4-dimensional gravity can be investigated. We are particularly interested 
on massive gravity on the Bloch brane, mainly on the effect of the massive spectrum of 
excitations for the localization of gravity on the brane. The form of the potential shows 
that besides the zero mode there are no other normalized states, with the massive states 
escaping from the brane as plane waves. 
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Figure 1: Plots of the Schroedinger-like potential Vsch{z) for r = 0.05 (thinner trace) and r = 0.10. 
Case r — 0.30 (thicker trace) corresponds to 1/2 Vsch{z) to ease comparison. 

The zero mode is responsible for gravity localization on the brane, whereas the massive 
modes introduce modifications on the Newtonian potential where in some scale gravity on 
the brane turns out to be essentially 5-dimensional. 

We numerically studied the Newtonian potential for two unit test masses separated by 
a distance R. In order to obtain the Newtonian potential one sums the tower of Kaluza- 
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Klein excitations to the usual contribution coming from the zero- mode as |13] 



1 1 r°° p-mR 
U{R) = G- + —J^ dm^|V^„(0)p, (2.7) 



where the 4-dimensional coupling G = M^^, is the fundamental 5-dimensional Planck 
scale and the integration is the considered at the brane position z = on the thin brane 
case. Separating from the Einstein-Hilbert action the 4-dimensional part from the extra 
dimensional one leads to an expression relating the two scales as |13| 



f + OO 

Ml = Ml / dze-3^(^)/2^ (2. 



In this way one can see that those scales are related to the integral of the squared wave- 
function for the zero- mode - as we have imposed ^/'o(O) = 1 we cannot associate a true 
probability for tI)q{z). That is the reason one can say integrability of the zero-mode is 
equivalent to a non-null 4-dimensional gravitational constant. Note that this does not 
guarantee localization of gravity in all scales, as the massive modes can modify sensibly 
the gravitational potential on the brane. However, see from Eq.(2.7) that the higher mas- 



sive modes are exponentially suppressed due to the Yukawa coupling. In this way gravity 
localization and the reproducibility of 4-dimensional gravity on the brane is strongly re- 
lated to the decoupling of the soft KK modes on the brane. In order to make easier the 
numerical approach to obtain the function U{R), see that the integral on the continuum 
modes is already saturated at m ~ 10/ R. 

Eq.(2.8) shows that |^m(0)P is of fundamental importance in order to achieve a true 



description of the Newtonian potential. As the massive modes are on the form of plane 
waves, normalization can be considered in a box of fixed length. The conformal transfor- 
mation from y to z is done in order to achieve a distribution of points on the z variable 
with fixed step with a controlled precision. With the function y{z) one can obtain A{z) and 
the Schroedinger potential Vgchiz)- The massive modes are obtained using the Numerov 
method |15], whose efficiency was evaluated after comparison of the m = case with the 



expression given by Eq. (2.6) for the zero-mode. For each value of the parameter r we 



have chosen a box extending from z — ~Zmax 

to •2- — Zmax ) with Zjnax chosen for each 
r such that the potential Vgch was near to the asymptotic value zero. This and a suffi- 
ciently large number of points turns out to be necessary in order to obtain better efficiency 
for comparison with the zero mode. We also expect that the normalization procedure for 
higher massive modes turns out to have higher confidence as the box collects larger number 
of wavelengths. We conclude from Fig. |l| that for lower values of r one needs boxes with 
larger z^ax ■ The Numerov method consists of a discretization of the Schroedinger equation 
in order to get a recurrence formula after choosing a fixed value for the non-normalized 
wavefunction at z = 0. 
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r 


Ml/M^ 


(Eq. ^ 


0.05 


14.67 


0.0681 


0.10 


8.85 


0.113 


0.15 


6.62 


0.151 


0.20 


5.43 


0.184 


0.30 


4.13 


0.241 


0.50 


3.04 


0.329 



Table I: Ratio between 4-dimensional coupling and fundamental 5-dimensional Planck scale 
for some values of r. Note the greater importance for Mj^ for smaller values of r. Third column 



contains the normalized zero-mode on the brane obtained using Eq. 2.6. This can be compared 
with the corresponding values for m — obtained with Numerov method and displayed on Fig. ^ 



Fig. shows the normalized probability for massive modes on the brane as a function of 
the masses of the modes for r = 0.3 and r = 0.05. First note that in the limit m — > 0, for the 
chosen parameters displayed on the same figure, Numerov method gives iV'oCO)!^ = 0.247 
for r = 0.3 and |^o(0)P = 0.0690 for r = 0.05. Comparing those values with Table I, this 
results in a relative error for determining the zero-mode of 2.5% and 1.3%, respectively. 
This shows that our method has higher precision for lower values of r. In particular, 
for higher values of r this error is difficult to reduce considerably as the Schroedinger 
potential goes slowly to zero in comparison to the volcano-like potential for lower values 
of r. Note also that the near- zero mode scattering states occur on the brane at higher 
probability, with the higher massive modes tending to a fixed probability. The asymptotic 
behaviour for higher massive modes shows that exponential suppression on Eq.( |2.7| ) is 
effective. The appearance of resonances could be identified as isolated peaks on energies 
below the maximum peak of the volcano potential and be responsible for a modification 
of the Newtonian potential obtained further in this section. From the figure one can see 
few evidence of resonances for r = 0.3. However, for r = 0.05 there is an oscillation on 
the curve for small masses, indicating the presence of resonances for values of r smaller 
than a critical one. In this way we see that the splitting of the Schroedinger potential for 
lower values of r, besides leading to an internal structure Q on the brane, also leads to 
a richer distribution of the massive modes. Extensions of Fig. |2| for higher masses up to 
m = 1000 show that the probabilities tend to constant values on both cases, indicating 
similar contributions to the 4-dimensional gravity on short distances, as we will see in the 
following. 

To better quantify the results, we are interested on the L parameter, defined as the 
exponent of the R-dependence of the Newtonian potential, 

U{R) = A, (2.9) 

where A is independent of and is related to the 5-dimensional gravitational constant G^. 
From the former equation we get L = — (9(log U) /9(log R). The dependence of the exponent 
L as a function of log R has an interesting interpretation. It describes the dimensionality 
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Figure 2: Normalized probability for finding massive modes on the brane as a function of the mass 
of the modes, for r = 0.3 (left) and r = 0.05 (right). Zmax is for box thickness, N the number of 
points and 6z the step on z variable. 



of the gravitational interaction as a function of the scale of separation between two test 
masses. 

Note, however, from Eq.( |2.7| ) that in order to obtain the Newtonian potential, an 
important information is the relation between the two scales M| and M^. Table I displays 
some results for their ratio depending on the r parameter. Note from the Table that the 
4-dimensional coupling constant G = I/M4 is more pronounced for larger values of r, 
indicating that the gravity localization is favored for those parameters. In particular, we 
see that when the field x is tuned to zero for the limit r = 0.5 we see that G achieves 
the maximum value compared to the fundamental scale. In this way we can associate the 
introduction of the additional field x the consequently richer structure of the brane to 
a larger modification of 4-dimensional gravity due to the extra dimension. 

Some results for L as a function of log R can be seen in Fig. ^ for coupling parameter 
r = 0.05 and r = 0.3. From the figure one can see that the exponent L interpolates between 
the usual 4-dimensional gravity L = 1 to a pure 5-dimensional gravity L = 2. The extra- 
dimensional characteristic is revealed on the brane for short separation between the test 
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masses, whereas the usual Newtonian potential is recovered for the whole scale of distances 
above a certain transition region. In Fig. ^ one can see that the explicit contribution for 
L coming from the continuum modes would induce a transition region with higher and 
lower values around the L = 2 region. In this way one can see that in the present model 
the normalized zero mode is essential for a proper recovering of the 4-dimensional gravity 
on the brane. Note also that on the transition region, one gets a higher dimensionality 
for lower values of r. If we compare this with Fig. |l| we see that the splitting of the 
Schroedinger potential favors the extra dimension to be revealed on the brane. 

Note that the Bloch brane scenario leads to different limit than RS ^ model for the 
Newtonian potential on short distance scales. As one knows, the AdS case for RS model 
gives L = 3 for those scales. However, the continuum modes decouple and the zero-mode 
is normalized in both scenarios, leading to corrections to Newton's law that are suppressed 
on large distance scales. 




Figure 3: Gravity localization on the brane for r ~ 0.3 and r = 0.05: for large scales the Newtonian 
potential is essentially 4-dimensional. 



3. Conclusions 

In this work we have studied the Bloch brane, with particular interest on signals from the 
extra dimension on the Newtonian potential. The spectrum analysis for the massive modes 
indicates a richer structure for low values for the parameter r, where some resonances do 
appear. This richer structure can be related to the splitting of the Schroedinger potential 
observed on Ref. Q. Also the resonance structure was revealed due to an improvement of 
the numerical normalization of the massive modes in a box and increased the understanding 
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of the Bloch brane. Our results for the ratio M^/M^ indicate that modifications coming 
from the extra dimension turns to be more important for lower values of the r parameter. 
This observation was confirmed from the Newtonian potential analysis where we have 
considered the change of the potential decay law U ~ with the increase of the 

separation R between two unit masses. We analyzed numerically the continuous change 
of the exponent L with log(i?) indicating the change of a pure 4-dimensional character for 
the gravitational interaction on the brane to a 5-dimensional one as one moves to shorter 
separation between the masses. In this way our result addresses the effort for experimental 
determination of gravitational interaction on short distances. In the Bloch brane scenario 
experimental modifications for the usual Newtonian law would indicate the presence of the 
extra dimension, for the parameter r very small. 
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